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Introduction 


Prostate  cancer  is  already  the  most  common  malignancy  in  men,  and  with  an  ever-aging  population 
the  incidence  of  the  disease  is  certain  to  increase.  Surgery  is  the  most  common  form  of  treatment  for 
localized  disease  but  many  thousands  of  patients  undergo  radiation  therapy  each  year  in  the  United 
States.  The  advent  of  PSA  testing  and  an  increasing  public  awareness  of  the  disease  suggest  that 
more  and  more  patients  will  select  radiation  treatment  to  avoid  the  potential  side  effects  of  surgery. 
The  use  of  radiation  therapy  is  not  definitive,  and  with  between  30  and  70%  of  patients  suffering  an 
increasing  PSA  5  years  after  treatment,  improved  methods  of  patient  follow-up  are  needed. 

Magnetic  resonance  imaging  (MRI)  offers  the  most  accurate  staging  information  when  cancer  is 
confined  to  the  gland.  However,  only  one  study  to  date  has  looked  at  patients  following  radiation 
treatment  and  the  findings  were  inconclusive.  A  preliminary  prospective,  quantitative  study,  such  as 
this,  is  an  essential  step  in  the  assessment  of  MRI  of  the  irradiated  gland.  The  information  obtained 
from  such  a  study  will  provide  important  physiological  and  structural  information  about  the  prostate 
gland  and  the  changes  induced  by  treatment.  These  changes  have  implications,  not  only  for  MR 
interpretation,  but  also  for  pathological  analysis  of  biopsy  specimens. 

The  following  report  details  the  final  results  of  a  quantitative  magnetic  resonance  imaging  (MRI)  study 
of  prostate  cancer  before  and  after  radiation  therapy.  The  study  aimed  to  examine  the  influence  of 
external  beam  radiation  on  the  prostate  and  prostate  cancer  using  novel  quantitative  MRI  techniques. 
Twenty-two  men,  previously  diagnosed  with  prostate  cancer,  were  studied  using  T1  and  12  magnetic 
resonance  relaxation  mapping  and  contrast  agent  kinetic  methods  before  and  after  treatment  by 
radiotherapy.  The  MRI  findings  will  be  correlated  with  biochemical  (prostate  specific  antigen)  and 
biopsy  results. 

Body 

The  following  section  summarises  the  research  accomplishments  associated  with  the  tasks  outlined  in 
the  approved  Statement  of  Work.  Despite  the  delays  experienced  in  transferring  the  grant  funding 
from  the  University  of  Florida  to  the  University  of  Manchester  and  the  subsequent  delay  in  obtaining 
Human  Subjects  approval  (detailed  in  the  Year  1  annual  report),  significant  progress  has  been  made 
with  the  Statement  of  Work.  The  research  accomplishments  associated  with  each  of  the  tasks 
outlined  are  described  below. 

Task  1 . 


Months  1-6.  To  complete  the  imaging  sequence  design  and  testing,  (specific  aims  1  and  2). 

a.  The  T1/T2  phantom  will  be  constructed. 

b.  Calibration  of  the  dynamic  imaging  sequences  (proton  density  weighted  +  T1 -weighted  sequences) 

c.  Optimization  of  the  FSE,  T2  imaging  sequence. 

d.  Selection  of  dynamic  and  T2  imaging  sequences.  Finalize  imaging  protocol. 


As  noted  in  the  Year  1  report,  with  the  cooperation  of  Drs.  Steve  Blackband,  Amanda  Barry  (Post- 
Doctoral  Research  Assistant  employed  during  the  first  year  of  the  study)  and  Geoff  Parker  (Research 
Fellow  working  on  related  studies),  all  elements  of  Task  1  were  successfully  completed. 

Task  2. 


Months  7-12.  Dynamic  data  simulation  and  model  development  (specific  aim  1).  Software  development 
(specific  aims  1  and  2). 

a.  Compare  simulated  data  with  patient  data  (both  existing  and  data  being  acquired  in  task  3). 

b.  Error  propagation  analysis  of  models  under  test. 

c.  Software  development.  T2  maps  and  dynamic  data  analysis. 
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a.  & 

b.  All  elements  of  Task  2  were  completed  (as  outlined  in  the  Year  2  and  Year  3  annual  reports)  but 
this  work  continues  to  provide  the  stimulus  for  further  related  developments.  In  particular,  the 
analysis  of  errors  in  the  techniques  employed  has  led  to  further  work  on  the  assessment  of 
dynamic  T1  measurement.  The  positive  findings  of  this  work  stress  the  importance  of  arterial 
input  function  assessment  for  minimising  measurement  errors.  Such  measurements  are  made  in 
this  study: 

i.  Appendix  1.  An  abstract  entitled,  WT1  estimation  using  variable  flip  angle  spoiled  gradient 
echo  for  dynamic  contrast-enhanced  MRI:  Arterial  input  measurement  improves  accuracy 
in  the  presence  of  B1  error”  was  accepted  for  presentation  at  the  12th  Annual  Meeting  of 
the  International  Society  for  Magnetic  Resonance  in  Medicine  (Kyoto,  May  2004).  The 
abstract  describes  errors  associated  with  contrast  agent  measurements  and  how  these 
may  be  minimised  by  data  analysis  employing  concurrent  measurement  of  the  arterial 
input  function. 

c.  Software  for  the  analysis  of  the  MR  data  is  in  place  and  has  been  used  in  the  subsequent  tasks. 
Examples  of  the  results  obtained  were  presented  in  the  Year  2  report  (Figures  1,  2  and  3). 

Task  3. 


Months  7-18.  Patient  recruitment  (21  patients)  and  examination  for  pretreatment  phase  of  specific  aim  3. 

a.  In  collaboration  with  Pathology,  examine  recent  biopsy  data  for  potential  patients. 

b.  Select  patients  (in  collaboration  with  Mr.  Clarke)  for  the  study.  Obtain  written,  informed  consent. 

c.  Acquire  data  from  the  1 .5  T  system.  Backup  (optical  disks)  and  transfer  data  to  laboratory. 

d.  Stain  and  reanalyze  biopsy  specimens. 


a.  & 

b.  A  total  of  22  patients  were  successfully  identified,  recruited  and  consented  by  the  end  of  February 
2003.  A  further  suitable  patient  was  included  in  the  study  after  writing  of  the  2003  annual  report. 

c.  Data  were  successfully  acquired  from  all  patients.  Ten  patients  (described  in  the  Year  2  report) 
were  scanned  between  July  2001  and  March  2002;  the  remaining  12  patients  were  scanned 
between  May  2002  and  March  2003.  All  data  were  backed-up  on  both  optical  disk  and  CD. 

d.  As  noted  in  the  Year  2  and  Year  3  reports,  it  had  been  decided  that  the  most  efficient  and 
appropriate  means  of  analysing  the  biopsy  specimens  was  to  process  them  as  a  batch.  However, 
it  has  since  become  apparent  that  we  may  receive  fewer  samples  than  originally  expected  thus 
work  is  now  underway,  by  Drs.  West  and  Clarke,  to  assess  the  possibility  of  analysing  pre¬ 
treatment  samples  obtained  from  various  outlying  hospitals  (the  original  referring  sites). 

Task  4 . 


Months  13-24.  Data  analysis,  pretreatment  phase. 

a.  Generation  of  Gd-DTPA  concentration  maps,  uptake  rate/max  uptake  maps  and  T2  maps. 

b.  Image  review  with  biopsy  results.  Preliminary  identification  of  cancer/BPH/normal  tissue. 

c.  Region  of  interest  analysis.  Model  fitting  (dynamic  data). 


a.  As  noted  in  the  Year  3  report,  maps  were  generated  for  each  patient.  Examples  of  these  were 
presented  in  the  Year  2  report  (Figs.  1 ,  2  and  3). 

b.  As  noted  in  the  Year  3  report,  Dr.  Charles  Hutchinson,  with  the  help  of  Mr.  Roberts  and  Dr. 

Buckley,  identified  regions  of  suspected  tumour,  normal  peripheral  zone  and  muscle  (internal 
obturator,  chosen  as  a  reference  tissue).  A  total  of  22  tumour  and  22  muscle  regions  were 
defined  in  22  patients.  Normal-appearing  peripheral  zone  was  not  visible  in  two  patients  thus  only 
20  normal  prostate  regions  were  defined. 

c.  As  noted  in  last  year’s  report,  these  data  were  all  analysed  successfully.  The  results  were 

subsequently  presented  at  the  11th  Annual  Meeting  of  the  International  Society  for  Magnetic 

Resonance  in  Medicine  (Toronto,  July  2003): 
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i.  Appendix  2.  An  abstract  entitled,  In  vivo  determination  of  the  microvascular 
characteristics  of  prostate  cancer  using  dynamic  contrast-enhanced  MRI”  was  presented, 
in  an  oral  clinical  science  focus  session  at  the  11th  Annual  Meeting  of  the  International 
Society  for  Magnetic  Resonance  in  Medicine  (Toronto,  July  2003).  The  abstract  (written  in 
November  2002)  details  only  the  first  few  patients  studied  but  the  entire  22-patient  data 
set  was  presented  at  the  meeting. 

These  data  have  now  been  fully  analysed  and  the  results  are  described  in  a  manuscript  that  has 
been  accepted  for  publication  by  the  journal  Radiology,  the  leading  publication  (highest  impact 
factor)  in  this  field  [1].  Full  details  of  these  results  are  presented  in  the  attached  manuscript 
(Appendix  3): 

ii.  Appendix  3.  A  full  paper,  “Vascular  characteristics  of  prostate  cancer:  Evaluation  with 
dynamic  contrast-enhanced  T1 -weighted  MRI  -  initial  experience”  in  press  in  the  journal 
Radiology  (accepted  for  publication  April  29th,  2004). 

In  addition  to  the  vascular  data,  the  relaxation  data  were  also  fully  analysed.  These  results 
(combined  with  the  complete  vascular  data  set)  are  shortly  to  be  presented  at  the  12th  Annual 
Meeting  of  the  International  Society  for  Magnetic  Resonance  in  Medicine  (Kyoto,  May  2004).  The 
results  are  detailed  in  the  abstract: 

i.  Appendix  4.  An  abstract  entitled,  “Characterization  of  prostate  cancer  by  quantitative  MRI” 
accepted  for  oral  presentation  at  the  12th  Annual  Meeting  of  the  International  Society  for 
Magnetic  Resonance  in  Medicine  (Kyoto,  May  2004).  The  abstract  describes  the 
complete  quantitative  data  set  from  the  pre-treatment  prostates. 

Task  5. 


Months  19-30.  Patient  follow-up  (21  patients)  and  examination  for  post-treatment  phase  of  specific  aim  3. 

a.  Request  2nd  examination  for  each  of  the  patients  studied  in  task  3. 

b.  Acquire  data  from  the  1.5  T  system.  Backup  (optical  disks)  and  transfer  data  to  laboratory. 


a.  & 

b.  To  date  19  patients  have  undergone  their  follow-up  MR  examinations  and,  as  previously,  all 
their  data  was  backed-up  to  both  optical  disk  and  CD.  The  mean  follow-up  time  was  14 
months  post-radiotherapy  (range  12  to  17  months).  Two  further  patients  are  due  follow-up  MR 
in  May  and  July  of  2004.  One  patient  developed  metastases  and,  due  to  their  illness,  was 
unable  to  undergo  follow-up  MR.  This  patient  has  since  died.  Dr.  Logue  and  his  staff  continue 
to  coordinate  the  follow-up  of  all  patients. 

Task  6 . 


Months  19-30.  Data  analysis,  post-treatment  phase. 

a.  Generation  of  Gd-DTPA  concentration  maps,  uptake  rate/max  uptake  maps  and  T2  maps. 

b.  Image  review  with  previous  MRI  and  biopsy  results.  Identification  of  cancer/BPH/normal  tissue. 

c.  Region  of  interest  analysis.  Model  fitting  (dynamic  data). 


a.  Relaxation  rate  maps  have  been  generated  for  19  patients  while  tracer  kinetics  analysis  is 
complete  on  12  patients. 

b.  The  task  of  identifying  cancer  and  normal  tissue  is  complicated  by  the  tissue  changes  caused  by 
radiotherapy.  These  include  generalised  loss  of  signal  from  the  prostate  gland  and,  in  particular, 
a  loss  of  internal  structural  contrast.  Dr.  Charles  Hutchinson  has  in  some  cases  identified  regions 
of  suspected  residual  tumour  while  in  others  has  simply  drawn  new  regions  to  match  as  closely 
as  possible  to  the  regions  drawn  on  the  pre-treatment  images.  A  similar  approach  has  been  taken 
for  the  normal  peripheral  zone.  Muscle  regions  have  also  been  identified. 

c.  Relaxation  data  from  19  patients  has  been  analysed  successfully.  The  results  are  outlined  in  the 
figure  below.  It  is  currently  too  early  to  speculate  on  the  significance  of  the  tracer  kinetics  analysis 
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(with  only  12  complete  data  sets  analysed  to  date)  and  a  full  analysis  will  be  postponed  until  the 
follow-up  study  is  complete. 


Figure  1. 

Relaxation  data  (T1  and  T2) 
measured  in  regions  of  tumour, 
normal  prostate  gland  and 
muscle.  Data  are  shown  for  19 
patients  examined  before  and 
14  months  after  radiotherapy. 
Significant  increases  following 
radiotherapy  are  seen  in 
tumour  T1  and  muscle  T2  (p  < 
0,01,  paired  t-test)  and  a  small 
decrease  is  seen  in  normal 
prostate  gland  T2  (p  <  0.05). 


Task  7. 


Months  31-36.  Review  of  results,  patient  follow-up  and  publication. 

a.  Obtain  follow-up  biopsies  on  each  patient  (18  months  after  treatment). 

b.  Stain  and  analyze  biopsy  specimens. 

c.  Review  each  patients  clinical  and  biochemical  data. 

d.  Compare  and  correlate  MRI  data  and  clinical  data. 

e.  Publication  of  significant  findings/experience. 


a.  To  date  follow-up  biopsies  have  been  requested  for  12  patients.  Three  patients  declined  to 
undergo  the  procedure  and  four  patients  have  relapsed  and  have  selected  further  treatment 
such  that  biopsy  was  not  recommended.  Thus  only  five  patients  have  undergone  follow-up 
biopsy. 

b.  Dr.  West  has  recently  taken  samples  from  the  five  follow-up  biopsies  and  is  assessing  them 
for  CD31/CD34  antibody  staining. 

c.  Dr.  Logue  and  his  team  continue  to  review  the  patients  studied.  We  have  complete  records  of 
all  biochemical  tests  (PSA)  and  the  patients  are  now  being  stratified  into  responders  and  non¬ 
responders.  We  anticipate  this  assessment  will  be  completed  by  the  summer. 

d.  The  comparison  of  histopathology,  clinical  findings  and  MR  data  is  now  underway.  The  pre¬ 
treatment  MR  data  will  be  compared  with  post-treatment  using  patient  response  (determined 
by  analysis  of  PSA  levels)  as  a  covariate.  These  data  will  be  reported  in  the  literature  as  soon 
as  they  have  been  fully  analysed  and  interpreted. 

e.  Two  full  papers  have  been  produced  from  the  data  accrued  and  analysed  thus  far  (one 
published  [2],  Appendix  5;  one  accepted  for  publication  [1],  Appendix  3).  It  is  anticipated  that 
at  least  two  further  papers  will  be  submitted  to  leading  journals  to  describe  the  findings  of  the 
post-treatment  data.  One  to  illustrate  the  changes  in  prostate  physiology  in  response  to 
radiotherapy  and  one  to  assess  the  prognostic  capabilities  of  contrast-enhanced  MRI  in  this 
patient  group. 

Key  Research  Accomplishments 


•  Implementation  of  a  quantitative  imaging  protocol  suitable  for  analysis  with  a  distributed 
parameter  tracer  kinetics  model  (i.e.  rapid  data  acquisition  with  arterial  input  measurement). 

•  Acquisition  of  data  from  22  patients  pre-treatment.  The  data  has  been  fully  analysed, 
assessed  and  published  in  a  leading  journal. 

•  Complete  follow-up  data  obtained  in  19  patients. 


7 


Final  Report  DAMD1 7-00-1-051 7 


•  Preliminary  analysis  of  relaxation  data  indicates  excellent  data  quality  and  post-radiotherapy 
changes  in  both  prostate  and  nearby  muscle. 


Reportable  Outcomes: 

The  outcomes  are  listed  in  chronological  order: 

1.  July  2003. 

Abstract.  G.J.  Parker,  A.  Jackson,  J.C.  Waterton,  D.L.  Buckley,  “Automated  arterial  input  function 
extraction  for  Tl-weighted  DCE-MRI”,  presented  at  the  11th  Annual  Meeting  of  the  International 
Society  for  Magnetic  Resonance  in  Medicine  (Toronto,  Canada). 

2.  July  2003. 

Abstract.  D.L.  Buckley,  C.  Roberts,  G.J.  Parker,  J.P.  Logue,  C.E.  Hutchinson,  “In  vivo  determination  of 
the  microvascular  characteristics  of  prostate  cancer  using  dynamic  contrast-enhanced  MRI”, 
presented  at  the  1 1th  Annual  Meeting  of  the  International  Society  for  Magnetic  Resonance  in  Medicine 
(Toronto,  Canada). 

3.  February  2004. 

Abstract.  D.L.  Buckley,  G.J.M.  Parker,  “T1  estimation  using  variable  flip  angle  spoiled  gradient  echo 
for  dynamic  contrast-enhanced  MRI:  Arterial  input  measurement  improves  accuracy  in  the  presence 
of  B1  error”  accepted  for  presentation  at  the  12th  Annual  Meeting  of  the  International  Society  for 
Magnetic  Resonance  in  Medicine  (Kyoto,  May  2004). 

4.  February  2004. 

Abstract.  D.L  Buckley,  C.  Roberts,  S.K.  Khaki,  G.J.M.  Parker,  J.P.  Logue,  C.E.  Hutchinson, 
“Characterization  of  prostate  cancer  by  quantitative  MRI”  accepted  for  oral  presentation  at  the  12th 
Annual  Meeting  of  the  International  Society  for  Magnetic  Resonance  in  Medicine  (Kyoto,  May  2004). 

5.  April  2004. 

Paper.  D.L.  Buckley,  C.  Roberts,  G.J.M.  Parker,  J.P.  Logue,  C.E.  Hutchinson  “Vascular 
characteristics  of  prostate  cancer:  Evaluation  with  dynamic  contrast-enhanced  Tl-weighted  MRI  - 
initial  experience”  accepted  for  publication  by  the  journal  Radiology. 


Conclusions 


The  study  is  nearing  completion  with  excellent  results  obtained  for  almost  all  aspects  of  the  work.  All 
patients  have  been  studied  successfully  prior  to  radiotherapy  and  these  data  form  the  basis  of  a 
manuscript  recently  accepted  for  publication  in  the  leading  journal,  Radiology  [1].  These  results  are 
both  novel  and  world  leading,  being  the  first  such  reported  measurements  using  MRI.  Recent  reports 
in  the  literature  highlight  the  developing  interest  in  this  field  of  MRI  [3,4]  and  another  CT  study  adopts 
a  similar  approach,  though  without  the  specificity  of  MRI  [5].  In  combination  with  the  vascular 
physiology,  MR  relaxation  measurements  match  well  to  previously  reported  results  [6].  All  that 
remains  to  be  assessed  in  the  pre-treatment  phase  is  the  biopsy  data,  and  the  effort  to  obtain  these 
data  is  on  going. 

The  follow-up  phase  of  the  study  is  almost  complete.  Only  one  patient  was  unable  to  complete  the 
study  having  developed  metastases.  Preliminary  analysis  of  the  relaxation  measurements  highlights 
the  quality  of  the  data  and  indicates  a  number  of  tissue  changes  14  months  after  radiotherapy.  We 
anticipate  further  insight  once  the  clinical  data  is  fully  interpreted  and  the  patients  are  stratified  into 
responders  and  non-responders.  At  this  point  the  prognostic  capabilities  of  the  MRI  techniques  will  be 
tested.  In  summary,  the  study,  despite  some  initial  delays,  has  been  a  resounding  success.  There 
remains  work  to  be  done  but  we  anticipate  future  studies  to  build  upon  these  successes. 
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Ti  estimation  using  variable  flip  angle  spoiled  gradient  echo  for  dynamic  contrast-enhanced  MRI:  Arterial  input  measurement 

improves  accuracy  in  the  presence  of  B,  error 


D.  L.  Buckley1,  G.  J.  Parker1 

imaging  Science  and  Biomedical  Engineering,  University  of  Manchester,  Manchester,  United  Kingdom 

Introduction. 

An  essential  step  in  the  quantitative  analysis  of  contrast  agent  uptake  is  the  assessment  of  changes  in  tissue  Ti.  The  rapid  acquisition 
times  often  required  for  these  measurements  dictate  that  they  are  made  using  gradient  echo  pulse  sequences  with  flip  angles  <  90°. 
Quantitative  measurements  made  using  such  sequences  rely  on  accurate  knowledge  of  Bi  and  it  is  well  known  that  many  transmit  coils 
can  produce  inhomogeneous  Bi  fields  and  this  can  be  compounded  by  slice  profile  and  dielectric  resonance  effects  [1,2].  In  this  study, 
the  effect  of  a  systematic  error  in  Bi  on  the  measurement  of  contrast  agent  uptake  in  tissues  is  examined. 

Methods. 


For  these  simulations  input  data  (1/Ti  against  time)  were  taken  from  regions  encompassing  a  prostate  tumour,  normal  peripheral  zone 
and  an  external  iliac  artery  (to  provide  an  arterial  input  function,  AIF)  following  the  administration  of  0.1  mmol/kg  Gd-DTPA-BMA.  The 
original  data  were  acquired  on  a  1.5  T  Philips  MR  system  every  2.3  s  for  4  minutes  using  a  3D  Ti-weighted  gradient  echo  pulse 
sequence  with  a  TR  of  2.5  ms  and  a  flip  angle  of  30°  [3].  Baseline  Ti  had  been  estimated  using  4  acquisitions  at  2°,  10°,  20°  and  30°  flip 
angles  [4].  From  the  temporal  changes  in  1/Ti  signal  intensity  time  courses  were  simulated  for  a  series  of  acquisitions  with  ±10%  and 
±30%  errors  in  Bi.  That  is,  for  acquisitions  with  true  flip  angles  of  21  °,  27°,  33°  and  39°  instead  of  the  nominal  30°.  These  errors  applied 
equally  to  the  data  used  for  baseline  Ti  estimates  (e.g.  a  -10%  error  in  the  flip  angle  will  result  in  acquisitions  at  1.8°,  9°,  18°  and  27° 
instead  of  the  nominal  2°,  10°,  20°  and  30°).  The  data  resulting  from  these  simulations  were  subsequently  converted  to  temporal 
changes  in  VT^  (using  the  erroneous  assumption  of  correct  flip  angles)  and  analyzed  using  a  two-compartment  model  providing 
estimates  of  Kirans ,  ve  and  vb  [5].  Two  analyses  were  performed:  1.  By  assuming  an  independently  measured  AIF  (errors  in  the  flip  angle 
do  not  affect  the  AIF  and  it  remains  in  its  original  form),  2.  By  assuming  the  AIF  is  measured  by  the  same  method  used  to  image  the 
tissue  (errors  in  the  flip  angle  affect  the  AIF  and  tissue  equally). 


Results. 


Fig.  1  True  1/Ti  time  course  in  an  artery  (faint  red  line,  scaled 
down  by  a  factor  of  10)  and  normal  prostate  (circles  with 
model  fit  in  red).  Estimated  1/Ti  time  course  in  artery  (faint 
blue  line)  and  prostate  (squares  with  model  fit  in  blue) 
when  flip  angle  error  =  -30%. 

The  true  baseline  Ti  of  the  prostate  and  artery  was  1000  ms  and 
1500  ms,  respectively.  A  -30%  error  in  the  flip  angle  results  in 
baseline  Ti  estimates  of  486  ms  and  728  ms,  respectively,  and  an 
overestimate  in  the  amplitude  of  the  1/Ti  time  courses.  Positive 
errors  in  the  flip  angles  led  to  overestimates  of  Ti  (up  to  1713  ms  for 
the  prostate)  and  concomitant  underestimates  in  the  amplitude  of 
the  1/Ti  time  courses.  When  the  simulated  data  were  presented  for 
fitting  by  method  1  errors  in  K'mns,  ve  and  vb  were  of  the  same 
magnitude  and  ranged  from  -24%  to  +40%  (resulting  from  +30%  to 
-30%  flip  angle  errors).  When  presented  for  fitting  by  method  2 
errors  in  K'rans  ranged  from  -3%  to  +10%;  errors  in  ve  ranged  from 
-3%  to  +1%  and  errors  in  vb  ranged  from  -6%  to  +21%. 


Discussion. 

Systematic  errors  in  the  amplitude  of  flip  angles  are  common  in  the  clinical  setting  [6j.  Correction  of  the  ensuing  errors  can  be 
accomplished  by  careful  system  calibration  [2,6]  but  this  may  be  impractical  in  many  situations.  When  encountered  in  dynamic  contrast- 
enhanced  studies  corrections  can  be  made  using  a  bookend  Ti  measurement  approach  [7].  However,  this  requires  accurate  baseline 
and  post-contrast  Ti  measurements  that  are  insensitive  to  Bi  errors.  When  left  uncorrected  errors  propagate  into  the  estimates  of  1/Ti 
(Fig.  1)  and  thereby  into  estimates  of  contrast  agent  concentration.  When  the  data  are  later  analysed  using  a  tracer  kinetic  model  the 
parameter  estimates  will  be  inaccurate.  Errors  in  K'rans  of  40%  result  from  a  -30%  error  in  Bi.  These  errors  are  magnified  as  TR 
increases  or  the  nominal  flip  angle  used  for  the  dynamic  acquisition  decreases.  However,  if  measurements  of  the  AIF  and  tissue  residue 
curve  are  made  simultaneously  then  subsequently  used  for  tracer  kinetic  analysis,  errors  in  the  resultant  parameter  estimates  are 
significantly  reduced.  Hence  a  -30%  error  in  Bi  results  in  only  a  10%  error  in  Ktrans.  This  improvement  arises  since  the  Bi  miscalculation 
affects  the  AIF  and  tissue  1/Ti  to  a  similar  degree.  Hence  the  use  of  a  measured  AIF  not  only  improves  the  underlying  accuracy  of 
parameter  estimates  [5]  but  also  reduces  sensitivity  to  Bi  error. 
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Synopsis. 

Tumour  angiogenesis  is  a  recognised  prognostic  indicator  in  patients  with  prostate  cancer  treated  with  external  beam  radiotherapy. 
Using  rapid  dynamic  contrast-enhanced  3D  MRI  and  a  distributed  parameter  tracer  kinetic  model  the  microvascular  characteristics  of 
prostate  cancer  were  assessed  in  a  group  of  patients  prior  to  treatment  with  radiotherapy.  Blood  flow  was  higher  in  the  tumours  than  in 
normal  peripheral  zone  but  there  were  no  significant  differences  in  PS  product  or  blood  volume.  This  novel  approach  may  prove 
valuable  as  a  prognostic  tool  and  the  patients  will  be  studied  post-treatment  to  assess  tissue  response. 


Introduction. 

Baseline  estimates  of  blood  flow  and  capillary  permeability  are  required  to  assess  response  of  the  prostate  to  treatment.  The  purpose  of 
this  study  was  to  assess  the  microvascular  characteristics  of  the  pathological  prostate  gland  prior  to  treatment  with  external  beam 
radiotherapy  using  dynamic  contrast-enhanced  MRI  and  a  distributed  parameter  tracer  kinetic  model  [1].  The  hypothesis  to  be  tested  is 
that  these  characteristics  will  be  sensitive  indicators  of  tissue  response  to  treatment  and  thereby  provide  a  useful  prognostic  tool. 


Methods. 

After  obtaining  written  informed  consent,  data  were  acquired  from  15  patients  with  biopsy  confirmed  stage  T2  or  T3  carcinoma  of  the 
prostate  gland.  These  patients  all  selected  external  beam  radiotherapy  for  subsequent  treatment  of  their  disease.  The  patients  were 
scanned  on  a  1.5  T  Philips  MR  system  using  a  pelvic  phased-array  coil  and  a  3D  Ti-weighted  gradient  echo  pulse  sequence.  Images 
were  acquired  using  flip  angles  of  2°,  10°,  20°  and  30°  to  estimate  baseline  Ti  [2].  This  was  followed  by  a  dynamic  series  in  which 
volumes  (flip  angle  30°)  were  acquired  every  2.3  s  for  approximately  4  minutes.  Early  in  this  series  0.1  mmol/kg  Gd-DTPA-BMA  was 
injected  at  3  ml/s  using  a  power  injector. 

The  image  data  were  subsequently  analysed  off-line.  For  each  patient  regions  of  interest  were  drawn  in  the  external  iliac  or  femoral 
arteries  (to  provide  a  vascular  input  function).  With  the  aid  of  T2-weighted  images  a  radiologist  (CEH)  drew  further  regions  in  prostate 
tumour,  muscle  (internal  obturator)  and,  where  possible,  normal  contralateral  peripheral  zone.  Signal  intensity  variations  in  these 
regions  were  converted  to  temporal  changes  in  Gd-DTPA-BMA  concentration  [2]  and  a  distributed  parameter  model  [1]  was  fitted  to  the 
data.  Parameters  estimates  in  each  region  were  compared  using  a  paired  t-test. 

Results. 

Analysis  of  the  data  from  each  region  provided  estimates  (mean  ±  SD)  of  blood  flow  (F),  blood  volume  (Vb)  and  microvascular 
permeability-surface  area  product  (PS)  [3]. 

Tissue  F  (ml/100  ml/min)  Vb  (ml/100  ml)  PS  (ml/100  ml/min) 

Tumour  (n=1 4)  81  ±66  2.1  ±3.7  24  ±11 

Normal  (n=12)  30  ±16  1.6  ±2.6  15  ±12 

Discussion. 

These  findings  highlight  the  capabilities  of  MR  imaging  to  estimate  parameters  only  previously  measured  using  PET  [4]  or  CT  [5].  The 
preliminary  results  confirm  that  blood  flow  to  tumour  tissue  exceeds  that  to  normal  prostatic  tissue  (p  =  0.02)  [4]  but  suggests  that  blood 
volume  and  capillary  permeability  may  not  be  significantly  different  (p  =  0.51  and  0.37,  respectively).  In  line  with  previous  MR  studies 
[6,7],  Kfrans  was  significantly  higher  in  tumour  than  normal  peripheral  zone  (p<0.0005).  Further  studies  are  ongoing  and  the  patients  will 
all  be  rescanned  12  months  after  treatment.  These  initial  findings  show  considerable  promise  for  isolating  the  physiological  basis  for 
contrast  enhancement  in  the  prostate. 


Figure  1.  Changes  in  1/T i  measured  overtime 
in  tumour  (left)  and  normal  peripheral  zone 
(right)  of  the  same  patient.  Blue  lines  show 
the  vascular  input  function  (scaled  down  by  a 
factor  of  10)  measured  in  the  external  iliac 
artery  and  the  red  lines  represent  model  fits  to 
the  data. 
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Abstract 


PURPOSE:  To  use  contrast-enhanced  magnetic  resonance  imaging  (MRI)  and  a  distributed 
parameter  tracer  kinetics  model  prospectively  for  evaluation  of  vascular  characteristics  of 
prostate  cancer. 

MATERIALS  AND  METHODS:  Twenty-two  patients  aged  between  57  and  76  years  old 
(mean,  67  years)  with  histologically  proven  adenocarcinoma  of  the  prostate  were  studied  at 
1.5  T  using  3D  dynamic  contrast-enhanced  ^-weighted  MRI.  Data  from  regions  of  interest 
drawn  in  tumor,  normal-appearing  peripheral  zone  tissue  and  muscle  were  analyzed  to 
provide  estimates  of  perfusion,  blood  volume,  interstitial  volume  and  microvascular 
permeability-surface  area  product.  These  estimates  were  compared  using  non-parametric 
Wilcoxon  Signed  Ranks  test. 

RESULTS:  Blood  flow  in  prostate  tumors  (n=22;  mean,  66  ml/100  ml/min)  was  significantly 
higher  than  in  contralateral  peripheral  zone  (n=20;  32  ml/100  ml/min).  Similarly  the  interstitial 
distribution  volume  in  tumors  was  enlarged  (42  ml/100  ml  in  tumor  versus  27  ml/100  ml  in 
normal  peripheral  zone).  Blood  volume  and  microvascular  permeability-surface  area  product 
in  tumors  (1.0  ml/100  ml  and  22  ml/100  ml/min,  respectively)  were  indistinguishable  from 
values  estimated  in  peripheral  zone  tissue  (1.5  ml/100  ml  and  21  ml/100  ml/min, 
respectively). 

CONCLUSION:  These  findings  show  considerable  promise  for  isolating  vascular 
characteristics  of  prostate  cancer. 

Keywords:  magnetic  resonance  (MR),  tissue  characterization;  prostate  neoplasms; 
neoplasms,  blood  supply. 


Introduction 


Magnetic  resonance  imaging  (MRI)  remains  the  most  promising  technique  for  the  detection 
and  staging  of  prostate  cancer  despite  its  current  limitations  in  sensitivity  and  specificity  [1]. 
Conventional  MRI,  typically  by  means  of  fast  T2-weighted  spin  echo  acquisitions,  provides 
excellent  anatomical  images  of  the  gland  in  which  the  normal  peripheral  zone  and  central 
gland  may  be  distinguished.  However,  malignant  tissue  is  not  always  clearly  visualized 
particularly  when  it  co-exists  with  benign  prostatic  hyperplasia  (BPH)  in  the  central  gland  [2]. 
Although  initially  considered  to  be  of  limited  use  for  studying  the  prostate  [3],  contrast- 
enhanced  MRI  has  shown  promise  in  a  series  of  studies  employing  rapidly  acquired  early  or 
dynamic  data  post-contrast  administration  [4-7].  In  these  studies,  cancers  (principally  in  the 
peripheral  zone  of  the  gland)  are  seen  to  enhance  more  rapidly  [5]  or  avidly  [7]  than  normal 
tissues. 

It  is  likely  that  the  early  enhancement  seen  in  cancer  is  due,  in  part,  to  angiogenesis 
associated  with  early  tumor  growth.  The  growth  of  prostate  cancer  is  associated  with  the 
development  of  a  rich  blood  supply  fed  by  a  large  network  of  immature,  leaky  blood  vessels 
[8],  The  density  of  this  network  is  related  to  tumor  grade  and  metastatic  potential  [9].  Since 
conventional  methods  of  assessing  angiogenesis  require  the  processing  of  pathological 
specimens  there  is  considerable  interest  in  developing  non-invasive  tools  to  assess  the 
microvascular  characteristics  of  prostate  cancer  in  vivo  [10].  A  link  between  contrast  uptake 
in  solid  tumors  and  their  microvascular  characteristics  certainly  exists  but  the  precise 
relationship  between  the  tracer  kinetics  measured,  often  in  terms  of  heuristic  measures,  and 
the  microvessel  density  measured  pathologically  remains  unclear  [11,  12].  The  assessment 
of  tracer  kinetics  in  absolute  physiological  terms  may  help  to  address  this  question  [13]. 
Increasing  microvessel  density  will  lead  to  an  increase  in  blood  flow,  blood  volume  and  the 
surface  area  of  vessel  walls.  An  upturn  in  vascular  endothelial  growth  factor  production  is 
likely  to  increase  the  permeability  of  these  vessel  walls.  Blood  flow,  blood  volume  and 
microvascular  permeability-surface  area  product  are  all,  in  principle,  quantifiable  through  the 
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analysis  of  contrast-enhanced  imaging  data  using  distributed  parameter  tracer  kinetics 
models  [10,  14].  These  have  yet  to  be  applied  in  vivo  to  a  meaningful  patient  group  using 
MRI,  though  a  recent  preliminary  X-ray  computerized  tomography  (CT)  study  has  shown 
promising  results  [15]. 

The  use  of  tracer  kinetics  modeling  in  the  analysis  of  tomographic  images  of  the  prostate  is 
widespread;  from  the  earliest  positron  emission  tomography  (PET)  study  [16],  through  a 
series  of  MRI  [17-19]  and  CT  studies  [15,  20].  While  the  MRI  studies  have  provided  sufficient 
spatial  and  contrast  resolution  to  identify  cancer,  normal  peripheral  zone  and  BPH  within 
individual  glands,  there  is  limited  information  about  the  microvascular  characteristics  of 
these  tissues.  Studies  employing  compartmental  models  have  shown  that  peripheral  zone 
cancers  have  higher  extraction-flow  products  (Kfrans,  [21])  than  normal  peripheral  zone  tissue 
[17,  18].  This  finding  is  supported  by  the  generally  higher  blood  flow  measured  in  tumor¬ 
bearing  prostates  compared  to  those  containing  BPH  alone  [16].  Differences  between 
tumors  and  normal,  BPH-containing,  tissue  in  the  central  gland  are  less  clear-cut  though  a 
recent  MRI  study  highlighted  potential  differences  in  flow-related  parameters  [22],  Thus,  the 
purpose  of  our  study  was  to  use  contrast-enhanced  MRI  and  a  distributed  parameter  tracer 
kinetics  model  prospectively  for  evaluation  of  vascular  characteristics  of  prostate  cancer. 
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Materials  and  Methods 


Tracer  kinetic  model 

There  are  distinct  advantages  to  employing  a  distributed  parameter  model  for  analysis  of 
tracer  kinetics  [23].  By  definition  compartmental  models  disregard  the  tracer’s  microvascular 
transit  time  and  it  is  assumed  to  be  negligible  on  the  time  scale  of  typical  Trweighted  MRI 
acquisitions.  Once  the  temporal  resolution  of  the  experiment  is  comparable  to  the  transit 
time,  the  use  of  a  distributed  parameter  model  becomes  a  realistic  option  [24].  Johnson  and 
Wilson  introduced  a  simple  distributed  parameter  model  in  which  the  internal  capillary  space 
contains  a  tracer,  the  concentration  of  which  varies  in  both  time  and  space  (along  the 
capillary  length)  [25].  Tracer  in  the  capillary  space  is  in  exchange  with  an  external  interstitial 
space  that  behaves  as  a  well-mixed  compartment  (i.e.  there  is  no  spatial  variation  in  tracer 
concentration).  St.  Lawrence  and  Lee  introduced  an  adiabatic  approximation  to  this  model  in 
1998  providing  a  time-domain  solution  suitable  for  the  analysis  of  imaging  data  [24]: 

C,is{t)  =  F  jCb{t-u)du  +  EFjch{u)exp  du  (1) 

0  x  V  ) 

where  Ctis  and  Cb  are  the  concentrations  of  contrast  agent  over  time,  t,  in  the  tissue  of 
interest  and  blood,  respectively.  F  is  the  tissue  blood  flow,  hct  is  the  blood  hematocrit  and  % 
is  the  mean  microvascular  transit  time  equal  to  the  blood  volume,  Vb,  divided  by  F.  Ve  is  the 
interstitial  tracer  distribution  volume,  u  is  a  dummy  integration  variable  and  the  extraction 
fraction,  E  =  1  -  exp(-PS/F(1-hct))  where  PS  is  the  microvascular  permeability-surface  area 
product.  The  product  EF(l-hct)  is  also  referred  to  as  K1™5  [21].  This  model  contains  4 
unknown  parameters  that  may  be  estimated  by  curve  fitting:  F,  E,  Ve  and  t.  From  these,  the 
parameters  PS  and  Vb  may  be  evaluated.  This  model  has  since  been  used  successfully  for 
the  analysis  of  both  MRI  and  CT  data  [15,  26,  27]. 

Patients 

Twenty-two  consecutive  men  with  histologically  proven  adenocarcinoma  of  the  prostate  who 
were  to  be  treated  with  neoadjuvant  hormonal  treatment  and  conformal  radiotherapy  were 


recruited  into  the  study  between  June  2001  and  February  2003.  The  Local  Research  Ethics 
Committee  approved  the  study  and  written  consent  was  obtained  from  all  men.  Histological 
diagnosis  was  obtained  by  transrectal  ultrasound  guided  biopsy  in  21  and  from  tissue 
obtained  at  the  time  of  a  transurethral  resection  of  the  prostate  in  one.  The  clinical  stage  at 
presentation  was  Tic  -  T3b,  the  modal  T-stage  being  T3  in  14  patients.  All  men  had  a 
negative  isotope  bone  scan.  The  mean  Gleason  score  was  6  (range,  5  -  8);  mean  PSA  27.5 
ng/ml  (range,  6-74  ng/ml);  and  mean  age  67  years  (range,  57  -  76  years).  Following  MRI 
and  in  accordance  with  local  protocol  all  men  were  treated  with  neoadjuvant  hormonal 
therapy  using  a  luteinizing  hormone  releasing  hormone  analogue  for  3-4  months  and  3D 
conformal  radiotherapy. 

Magnetic  resonance  imaging 

MRI  was  performed  at  1.5  T  (Gyroscan  NT/Intera;  Philips  Medical  Systems,  Best,  the 
Netherlands)  using  a  4-element  phased-array  receiver  coil  wrapped  around  the  pelvis. 
These  men  were  expected  to  undergo  follow-up  MRI  after  radiotherapy  thus  use  of  an 
endorectal  receiver  was  precluded  [28].  Following  the  acquisition  of  scout  images  data  were 
acquired  encompassing  the  abdomen  and  pelvis  in  both  the  axial  and  coronal  planes  using 
^-weighted  spin-echo  sequences.  These  images  were  used  to  assess  the  possibility  of 
extra-prostatic  spread.  Subsequently,  T2-weighted  fast  spin  echo  images  of  the  prostate 
were  obtained  in  the  sagittal  and  axial  planes  (field  of  view  =  200  x  200  mm,  slice  thickness 
=  3  mm,  gap  =  0.3  mm,  matrix  =  256  x  190,  repetition  time  (TR)  =  4111  ms,  effective  echo 
time  (TE)  =  135  ms,  echo  train  length  19  and  6  averages).  Following  the  acquisition  of 
quantitative  data  during  contrast  agent  administration  (detailed  below),  post-contrast  Tr 
weighted  spin  echo  images  of  the  abdomen  and  pelvis  were  again  acquired  in  the  axial 
plane. 

Quantitative  data  encompassing  the  pelvic  region  from  the  apex  of  the  prostate  gland  into 
the  bladder  were  acquired  using  an  axial  3D  spoiled  gradient  echo  sequence  0VFFE  with 
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an  elliptical  k-space  shutter;  field  of  view  =  300  x  300  x  100  mm,  matrix  following 
interpolation  =  128  x  128  x  20,  TR  =  2.5  ms,  TE  =  0.86  ms).  Prior  to  contrast  agent 
administration  the  native  Ti  relaxation  time  was  mapped  using  a  series  of  volume 
acquisitions  with  flip  angles  of  2,  10,  20  and  30  degrees  (each  acquired  with  4  averages). 
Subsequently,  a  dynamic  series  was  acquired  using  a  30-degree  flip  angle  and  no  averaging 
over  a  period  of  ~4  minutes  with  a  total  of  100  to  106  volumes  acquired;  one  every  2.3  s. 
Contrast  agent  was  administered  to  the  patients  through  a  cannula  placed  in  the  antecubital 
vein  prior  to  scanning.  Following  the  acquisition  of  the  first  5  volumes  of  the  dynamic  series 
0.1  mmol/kg  Gd-DTPA-BMA  (Omniscan;  Amersham-Health,  Little  Chalfont, 
Buckinghamshire,  UK)  was  injected  at  3  ml/s  using  an  MRI-compatible  power  injector 
(Spectris;  Medrad,  Indianola,  PA).  This  was  immediately  followed  by  a  similar  volume  flush 
of  normal  saline. 

Data  collection  and  statistical  analysis 

Data  were  transferred  offline  for  quantitative  analysis.  For  each  patient  4-pixel  (22  mm2) 
regions  of  interest  were  placed  (by  CR)  in  the  external  iliac  or  femoral  arteries  (to  provide  a 
vascular  input  function).  With  the  aid  of  the  biopsy  results  and  T2-weighted  images  a 
radiologist  (CEH)  drew  further  regions  in  muscle  (internal  obturator,  selected  as  a  control 
tissue;  region  size  range  148  to  956  mm2),  hypointense  areas  within  the  peripheral  zone  that 
were  suspected  of  containing  prostate  tumor  and,  where  possible,  normal  contralateral 
peripheral  zone.  These  regions  were  drawn  to  encompass  as  much  tumor  and  normal 
contralateral  peripheral  zone  as  was  visible  in  the  selected  slice.  Signal  intensity  variations 
in  these  regions  were  converted  to  temporal  changes  in  Gd-DTPA-BMA  concentration  using 
estimates  of  native  ^  calculated  using  the  pre-contrast  data  [29,  30].  Finally,  the  distributed 
parameter  model  described  by  St.  Lawrence  and  Lee  [24]  was  fitted  to  the  data.  Multiple  fits 
were  performed  to  each  tissue  region  using  a  range  of  starting  points  (systematic  variation  in 
the  initial  estimates  of  -c).  The  final  solution  selected  was  that  which  minimized  the  sum  of 
squared  differences  between  the  data  and  the  model  fit  [26,  27].  To  determine  if  differences 
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existed  between  vascular  characteristics  of  tumor,  normal  prostate  and  muscle  in  the 
patients  studied  parameter  estimates  were  compared  using  the  non-parametric  Wilcoxon 
Signed  Ranks  test  (SPSS  for  Windows,  version  1 1 .5,  2002).  Statistical  significance  was 
judged  at  p  <  0.05. 
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Results 


Magnetic  resonance  imaging 

Data  were  successfully  acquired  from  all  22  patients.  Hypointense  lesions  were  outlined  in 
each  patient  with  reference  to  the  T2-weighted  images  and  confirmed  by  reference  to  the 
biopsy  reports.  Normal  appearing,  hyperintense  regions  were  further  outlined  in  20 
prostates.  Example  data  are  presented  in  Figs.  1  to  3. 

Data  Analysis 

Analysis  of  the  data  from  the  tissue  and  arterial  regions  provided  estimates  of  the  relaxation 
and  tracer  kinetics  parameters.  Their  mean  values  ±  SD  are  quoted  in  the  text  below.  The 
native  Ti  relaxation  times  of  tumor,  normal-appearing  peripheral  zone  and  muscle  were  916 
±  298  ms,  962  ±  273  ms  and  1026  ±  242  ms,  respectively.  The  model  parameter  estimates 
did  not  follow  a  normal  distribution  and  the  range  and  median  values  for  each  region  and 
parameter  are  quoted  in  Table  1.  Blood  flow  in  the  tumors,  66  ±  43  ml/min/100  ml,  was 
significantly  greater  than  blood  flow  measured  in  normal-appearing  peripheral  zone,  32  ±  36 
ml/min/100  ml,  (p  <  0.001)  and  muscle,  9  ±  7  ml/min/100  ml,  (p  <  0.001).  Furthermore,  the 
interstitial  volume  of  contrast  agent  distribution  in  tumors,  42  ±  20  ml/100  ml,  was  larger  than 
that  measured  in  both  normal-appearing  peripheral  zone,  27  ±  10  ml/100  ml,  (p  =  0.003)  and 
muscle,  12  ±  17  ml/100  ml,  (p  =  0.002).  The  permeability-surface  area  product  of  vessels  in 
the  tumors,  22  ±  12  ml/min/100  ml,  was  similar  to  that  measured  in  normal-appearing  tissue, 
21  ±  24  ml/min/100  ml,  but  was  higher  than  seen  in  muscle,  5  ±  5  ml/min/100  ml,  (p  <  0.001). 
These  trends  were  lost  for  estimates  of  the  blood  volume  where  significant  differences  were 
not  seen.  Blood  volume  estimates  in  tumors,  1.0  ±  1.4  ml/100  ml,  normal-appearing 
peripheral  zone,  1.5  ±  2.4  ml/100  ml,  and  muscle,  1.8  ±  2.0  ml/100  ml,  were  all  very  low.  The 
combined  parameter  extraction-flow  product  (K,ran7(1-hct))  estimated  in  tumors,  28  ±  13 
ml/min/100  ml,  was  significantly  greater  than  in  normal-appearing  peripheral  zone,  15  ±  8 
ml/min/100  ml,  (p  <  0.001)  and  muscle,  5  ±  3  ml/min/100  ml,  (p  <  0.001). 
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Discussion 


These  findings  demonstrate  that  MRI  can  be  used  to  estimate  parameters  only  previously 
measured  in  the  prostate  gland  using  PET  [16]  or  CT  [15].  The  results  confirm  that  blood 
flow  to  tumor  tissue  exceeds  that  to  normal  prostatic  tissue  but  that  the  difference  in  blood 
volume  is  insignificant  [16].  An  increase  in  the  interstitial  distribution  volume,  Ve,  was 
measured  in  tumors  [18].  This  may  explain  the  high  peak  enhancement  measured  in  tumors 
[7]  and  reflect  disruption  of  the  epithelial  cell  layer  surrounding  prostatic  ducts  caused  by 
invading  cancer  cells  [31].  Unlike  previous  quantitative  MRI  studies  3D  data  were  acquired 
encompassing  the  entire  prostate  gland.  This  had  the  additional  advantage  of  enabling 
measurement  of  a  vascular  input  function  close  to  the  tissue  [30]. 

Validation.  Our  estimates  of  native  tissue  Ti  made  using  variable  flip  angle  3D  gradient  echo 
acquisitions  [29]  compare  well  with  previous  measures  in  both  muscle  [32]  and  the  prostate 
[33].  Single  capillary  distributed  parameter  models  have  previously  been  employed  for  the 
analysis  of  PET  data  [23]  and  the  St.  Lawrence  and  Lee  model  specifically  has  been  used  in 
both  MRI  [26]  and  CT  [15]  studies.  A  variant  of  the  model  has  been  validated  in  two  separate 
studies  for  the  measurement  of  blood  flow  using  CT  by  comparison  with  microsphere 
measurements  [34,  35].  Though  not  formally  validated  in  the  prostate,  simulation  studies 
have  indicated  that  the  St.  Lawrence  and  Lee  model  performs  well  with  comparison  to  both 
an  independent  multiple  path  distributed  parameter  model  and  commonly  used 
compartmental  models  [27],  In  line  with  previous  work  [17,  18],  the  combination  of  E  and  F 
(K"ans)  measurec|  jn  the  current  study  was  found  to  be  significantly  higher  in  tumor  than 
normal  peripheral  zone.  The  estimates  of  blood  flow  compare  well  with  estimates  made 
using  PET  [16].  Even  the  small  decrease  in  blood  volume  that  was  associated  with  cancer 
follows  the  trend  observed  in  that  earlier  study  [16],  Finally,  the  model-based  estimates  of 
blood  flow  and  microvessel  PS-product  correlate  well  with  the  predictions  made  by 
pathologists.  While  tumors  contain  a  greater  number  of  vessels  than  normal  tissue, 
specifically  in  angiogenic  hot-spots  [8],  the  vessel  surface  area  density  in  the  tumor  is 
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actually  reduced  [36].  This  reflects  the  size  and  shape  of  tumor  microvessels  that  tend  to  be 
broad  and  less  branched  than  those  in  normal  tissue  or  low  grade  tumors.  This  reduced 
vessel  surface  area  in  the  high  grade  cancers  is  offset  by  increased  blood  flow  and  vessel 
wall  permeability  [36],  physiological  changes  that  are  obscured  in  the  estimate  of  PS-product 
and  may  explain  the  similarities  between  the  estimates  made  in  tumors  and  normal  tissue. 

Study  limitations.  Though  concordant  imaging  and  biopsy  results  provide  a  highly  specific 
marker  of  cancer  in  the  peripheral  zone  [37]  the  study  lacks  a  reference  standard  such  as 
whole-mount  prostatectomy  specimens  [38],  Essential  to  the  analysis  of  any  experiment  is 
data  quality.  The  data  acquired  from  the  external  iliac/femoral  arteries  and  the  prostate  gland 
were  of  good  quality  but  the  signal  obtained  from  the  internal  obturator  muscle  suffered  from 
a  poor  signal-to-noise  ratio  (SNR)  and  sporadic  motion  artifact.  Motion  interfered  with  a 
number  of  the  signal-time  courses  analyzed  and  future  studies  will  benefit  from  image 
registration  prior  to  data  analysis.  Conventionally  the  SNR  is  addressed  at  the  expense  of 
temporal  resolution.  This  was  not  an  option  in  this  study  since  it  was  necessary  to  acquire 
data  very  rapidly  so  that  a  distributed  parameter  model  could  be  used.  Indeed,  the 
acquisition  time  employed  is  already  longer  than  recommended  for  similar  studies  [39].  The 
limited  SNR  is  a  common  problem  in  prostate  MRI  and  may  be  partly  addressed  by  the 
application  of  an  endorectal  coil  [38].  Though  certain  to  improve  the  signal  obtained  from  the 
prostate  it  is  unlikely  to  increase  the  signal  from  the  artery  or  muscle.  Moreover,  the 
endorectal  coil  may  be  contraindicated  in  the  assessment  of  patients  following  radiotherapy, 
such  as  the  subjects  studied  here  [28].  Alternatively,  SNR  could  be  improved  by  the  use  of  a 
high-field  MRI  system  [40].  These  limitations  may  have  influenced  the  accuracy  and 
precision  of  the  parameter  estimates  obtained.  An  earlier  study  revealed  that  estimates  of 
blood  volume  obtained  using  this  implementation  of  the  St.  Lawrence  and  Lee  model  were 
low  [27]  and  the  current  results  fall  below  previous  estimates  made  using  PET  [16].  This  may 
reflect  noise  or  an  inadequate  temporal  resolution  in  the  experimental  data  [39]  that  leads  to 
limitations  in  definition  and  integration  of  the  arterial  input  [20],  It  may  also  indicate  the 
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existence  of  a  limited  transendothelial  water  exchange  rate  [41,  42],  The  influence  of  such 
effects  may  be  reduced  by  further  decreasing  the  sensitivity  of  the  imaging  sequence  to 
water  exchange  (by  increasing  the  flip  angle,  so-called  exchange-minimization  [41]).  Such 
an  approach  further  compromises  the  SNR  of  the  sequence  and  a  practical  trade-off  must 
be  achieved. 

Clinical  implications.  These  results  provide  new  quantitative  physiological  data  concerning 
the  prostate  gland.  The  methods  represent  an  advance  on  previous  MRI  studies  [7,  18]  with 
measurement  of  an  arterial  input  function,  3D  coverage  of  the  entire  prostate  gland  and 
absolute  estimates  of  both  flow  and  PS-product.  Moreover,  similar  CT  studies  lack  tissue 
coverage  (only  one  or  two  sections  may  be  examined)  and  the  intrinsic  tissue  contrast  of 
MRI  that  enables  discrimination  of  tumor  and  normal  peripheral  zone  tissue  [15,  20].  The 
techniques  will  allow  further  studies  to  characterize  tissue  types,  improve  diagnosis  and 
categorize  patient  prognosis.  The  promising  results  obtained  using  qualitative  or  heuristic 
methods  could  be  assessed  to  identify  their  underlying  physiological  basis  [7,  22],  Tissue 
characterization  in  the  current  study  was  limited  by  the  method  of  region  selection  dictated 
by  the  patient  group  examined.  Nevertheless,  concordant  imaging  and  biopsy  results 
providing  a  highly  specific  marker  of  cancer  in  the  peripheral  zone  [37]  guided  selection  of 
prostate  tumor  regions.  Future  work  would  benefit  from  comparison  with  whole-mount 
prostatectomy  specimens  providing  an  improved  reference  standard  [38]. 
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Figure  captions 

Fig.  1  Tissue  enhancement  following  bolus  administration  of  Gd-DTPA-BMA.  Axial  T2- 
weighted  turbo  spin  echo  image  (TR/TE  4111/135  ms)  of  an  example  patient 
showing  the  prostate  gland  (a).  Regions  of  interest  have  been  drawn  to  outline  tumor 
(left/anterior  of  the  prostate)  and  contralateral  normal-appearing  peripheral  zone 
(right  side  of  the  prostate).  The  graphs  illustrate  enhancement  over  time  in  the 
external  iliac  artery  (faint  lines,  scaled  down  by  a  factor  of  10),  tumor  (b)  and  normal 
(c)  regions.  In  each  graph  the  bold  line  represents  the  best  fit  of  the  model  to  the 
data  and  provides  estimates  of  (amongst  other  parameters)  blood  flow,  F,  and 
interstitial  distribution  volume,  Ve.  (b)  Tumor  tissue.  F  =  109  ml/100  ml/min  and  Ve  = 
30  ml/100  ml.  (c)  Normal  tissue.  F  =  21  ml/100  ml/min  and  Ve  =  34  ml/100  ml. 

Fig.  2  (a)  Axial  T2-weighted  turbo  spin  echo  image  (TR/TE  4111/135  ms)  of  a  second 
example  patient  showing  the  prostate  gland  with  tumor  in  the  left  posterior  area  of  the 
gland,  (b)  Graph  illustrates  enhancement  over  time  in  tumor  (open  squares),  normal 
(filled  circles)  and  muscle  (open  circles)  regions.  The  bold  lines  represent  the  best  fit 
of  the  model  to  the  data.  Parameter  estimates  were  obtained  for  tumor,  F  =  64 
ml/100  ml/min  and  Ve  =  50  ml/100  ml;  normal,  F  =  19  ml/100  ml/min  and  Ve  =  40 
ml/100  ml;  muscle,  F  =  12  ml/100  ml/min  and  Ve  =  7  ml/100  ml. 

Fig.  3  (a)  Axial  T2-weighted  turbo  spin  echo  image  (TR/TE  41 1 1/135  ms)  of  a  third  example 
patient  showing  the  prostate  gland  with  tumor  in  the  left  posterior  area  of  the  gland, 
(b)  Graph  illustrates  enhancement  over  time  in  tumor  (open  squares),  normal  (filled 
circles)  and  muscle  (open  circles)  regions.  The  bold  lines  represent  the  best  fit  of  the 
model  to  the  data.  Note  the  noise  in  the  muscle  data  and  the  undulating  nature  of  the 
model  fits  reflecting  variability  in  this  patient’s  arterial  input  function.  Parameter 
estimates  were  obtained  for  tumor,  F  =  132  ml/100  ml/min  and  Ve  =  55  ml/100  ml; 
normal,  F  =  109  ml/100  ml/min  and  Ve  =  40  ml/100  ml;  muscle,  F  =  8  ml/100  ml/min 
and  Ve  =  18  ml/100  ml. 
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Introduction. 

Quantitative  estimates  of  the  baseline  physiology  and  MR  characteristics  of  the  prostate  may  help  to  predict  the  response  of  the  gland 
to  treatment.  The  purpose  of  this  study  was  to  assess  these  characteristics  in  the  pathological  prostate  gland  prior  to  treatment  with 
external  beam  radiotherapy  using  T-i  and  T2  mapping,  dynamic  contrast-enhanced  MR!  and  a  distributed  parameter  tracer  kinetic  model. 
The  hypothesis  to  be  tested  is  that  these  characteristics  will  serve  as  sensitive  indicators  of  tissue  response  to  treatment  and  thereby 
provide  a  useful  prognostic  tool. 

Methods. 

Twenty-two  men  with  histologically  proven  adenocarcinoma  of  the  prostate  were  recruited  into  the  study.  The  Local  Research  Ethics 
Committee  approved  the  study  and  written  consent  was  obtained  from  all  men.  Histological  diagnosis  was  obtained  by  transrectal 
ultrasound  guided  biopsy  in  21  and  from  tissue  obtained  at  the  time  of  a  transurethral  resection  of  the  prostate  in  one.  The  clinical  stage 
at  presentation  was  Tic  -  T3b.  All  men  had  a  negative  isotope  bone  scan.  The  mean  Gleason  score  was  6  (range,  5  -  8);  mean  PSA 
27.5  ng/ml  (range,  6  -  74  ng/ml)  and  mean  age  67  years  (range,  57  -  76  years).  These  patients  all  selected  external  beam  radiotherapy 
for  subsequent  treatment  of  their  disease.  The  patients  were  scanned  on  a  1.5  T  Philips  MR  system  using  a  pelvic  phased-array  coil. 
T2-weighted  fast  spin  echo  images  of  the  entire  gland  were  acquired  using  a  TR  of  4.75  s  at  echo  times  of  7,  45,  100  and  240  ms  to 
estimate  baseline  T2  [1].  A  3D  ^-weighted  gradient  echo  pulse  sequence  was  used  at  flip  angles  of  2°,  10°,  20°  and  30°  to  estimate 
baseline  Ti  [2].  This  was  followed  by  a  dynamic  series  in  which  volumes  (flip  angle  30°)  were  acquired  every  2.3  s  for  approximately  4 
minutes.  Early  in  this  series  0.1  mmol/kg  Gd-DTPA-BMA  was  injected  at  3  ml/s  using  a  power  injector. 

The  image  data  were  subsequently  analysed  off-line.  For  each  patient  regions  of  interest  were  drawn  in  the  external  iliac  or  femoral 
arteries  (to  provide  a  vascular  input  function).  With  the  aid  of  T2-weighted  images  a  radiologist  (CEH)  drew  further  regions  in  prostate 
tumour,  muscle  (internal  obturator)  and,  where  possible,  normal  contralateral  peripheral  zone.  For  each  region  estimates  of  Ti  and  T2 
were  made  and  signal  intensity  variations  were  converted  to  temporal  changes  in  Gd-DTPA-BMA  concentration  [2]  and  a  distributed 
parameter  model  [3]  was  fitted  to  the  data. 

Results. 

Analysis  of  the  data  from  each  region  provided  estimates  (mean  ±  SD)  of  Ti,  T2,  blood  flow  (F),  blood  volume  (Vb),  microvascular 
permeability-surface  area  product  (PS)  and  interstitial  volume  (Ve)  [4]. 

Ti  T2  F  Vb  PS  Ve 


Tissue _ (ms) _ (ms) _ (ml/1 00  ml/min)  (ml/1 00  ml)  (ml/1 00  ml/min)  (ml/1 00  ml) 


Tumour  (n=22) 

916  ±298 

93  ±15 

66  ±43 

1.0  ±1.4 

22  ±12 

42  ±20 

Normal  (n=20) 

962  ±  273 

131  ±27* 

32  ±  36* 

1.5  ±2.4 

21  ±24 

27  ±10* 

Muscle  (n=22) 

1026  ±  242 

52  ±8* 

9  ±7* 

1.8  ±2.0 

5  ±5* 

12  ±17* 

*Significantly  different  from  tumour  values  using  the  non-parametric  Wilcoxon  Signed  Ranks  test  (p<0.01). 


Discussion. 

Our  estimates  of  the  relaxation  times  compare  well  with  previous  results  and  underscore  the  drop  in  T2  seen  in  cancer  [1].  This  analysis 
was  based  upon  estimates  of  T2  made  using  only  the  first  three  echoes.  There  was  an  indication  of  the  multiexponential  T2  behaviour 
reported  by  Kjaer  et  al.  [5]  when  the  240  ms  echo  data  were  considered  and  this  finding  requires  further  study. 

The  results  of  the  tracer  kinetics  analysis  confirm  that  blood  flow  to  tumour  tissue  exceeds  that  to  normal  prostatic  tissue  but  that  the 
difference  in  blood  volume  is  insignificant  [6].  Similarly,  the  interstitial  distribution  space  in  the  tumours  is  enlarged  [7]  but  no  difference 
between  the  permeability-surface  area  products  of  tumour  and  normal  microvessels  was  detectable.  In  a  recent  publication  reporting 
the  use  of  contrast-enhanced  CT  and  the  same  tracer  kinetics  model  the  authors  stressed  the  importance  of  interpreting  results  with 
caution  since  the  precision  of  parameter  estimates  can  be  low  [8].  Indeed,  the  estimates  of  PS  made  in  this  study  were  not  always 
reliable;  PS  cannot  be  estimated  when  the  contrast  agent  is  completely  extracted  on  its  first  pass  (E  =  1).  The  blood  volume  estimates 
reported  by  Henderson  et  al.  [8]  were  larger  than  found  in  this  study  and  may  reflect  MR-specific  problems  associated  with 
transendothelial  water  exchange  [9]  in  addition  to  the  difficulties  previously  experienced  in  fitting  complex  models  [4].  However,  MR 
studies  such  as  this  have  distinct  advantages  over  CT  studies.  The  exquisite  soft  tissue  contrast  of  T2-weighted  images  allows  for  the 
analysis  of  hypointense  regions  against  normal-appearing  peripheral  zone  in  any  part  of  the  gland.  The  radiation  burden  of  the  dynamic 
CT  study  prevents  the  examination  of  more  than  one  or  two  sections  and  these  must  be  chosen  with  little  prior  knowledge  of  the 
location  of  cancer.  MR  may  be  used  repeatedly;  further  studies  are  ongoing  with  patients  returning  for  follow-up  MR  scans  12  months 
after  treatment.  These  findings  show  considerable  promise  for  isolating  both  the  MR  and  physiological  characteristics  of  prostate  cancer 
and  may  help  to  assess  tissue  response  to  treatment  in  this  and  future  studies. 

Acknowledgments.  Funded  by  the  U.S.  Department  of  Defense  Prostate  Cancer  Research  Program  (PC991154).  Jeanette  Lyons, 
Yvonne  Watson  and  David  Clark  kindly  helped  with  the  patient  studies. 

References.  1.  Liney  GP  et  al.  Magn  Reson  Imaging  15:1 177-1 186  (1 997);  2.  Zhu  XP  et  al.  J  Magn  Reson  Imaging  1 1 :575-585  (2000); 
3.  St.  Lawrence  KS,  Lee  TV  J  Cereb  Blood  Flow  Metab  18:1365-1377  (1998);  4.  Buckley  DL  Magn  Reson  Med  47:601-606  (2002); 

5.  Kjaer  L  et  al.  Magn  Reson  Imaging  5:23-30  (1987);  6.  Inaba  T  J  Urol  148:1457-1460  (1992);  7.  Padhani  AR  et  al.  Clin  Radiol  55:99- 
109  (2000);  8.  Henderson  E  etal.  Phys  Med  Biol  48:3085-31 00  (2003);  9.  Buckley  DL  Proc.  1Cfh  meeting  ISMRM,  Honolulu,  2120 
(2002). 


NOTES 


Magnetic  Resonance  in  Medicine  47:601-606  (2002) 
DOI  10.1002/mrm.l0080 


Uncertainty  in  the  Analysis  of  Tracer  Kinetics  Using 
Dynamic  Contrast-Enhanced  7^ -Weighted  MRI 

David  L.  Buckley* 


In  recent  years  a  number  of  physiological  models  have  gained 
prominence  in  the  analysis  of  dynamic  contrast-enhanced  T.,- 
weighted  MRI  data.  However,  there  remains  little  evidence  to 
support  their  use  in  estimating  the  absolute  values  of  tissue 
physiological  parameters  such  as  perfusion,  capillary  perme¬ 
ability,  and  blood  volume.  In  an  attempt  to  address  this  issue, 
data  were  simulated  using  a  distributed  pathway  model  of 
tracer  kinetics,  and  three  published  models  were  fitted  to  the 
resultant  concentration-time  curves.  Parameter  estimates  ob¬ 
tained  from  these  fits  were  compared  with  the  parameters  used 
for  the  simulations.  The  results  indicate  that  the  use  of  com¬ 
monly  accepted  models  leads  to  systematic  overestimation  of 
the  transfer  constant,  K^8"8,  and  potentially  large  underesti¬ 
mates  of  the  blood  plasma  volume  fraction,  Vp .  In  summary, 
proposals  for  a  practical  approach  to  physiological  modeling 
using  MRI  data  are  outlined.  Magn  Reson  Med  47:601-606, 2002. 
©  2002  Wiley-Liss,  Inc. 

Key  words:  perfusion;  gadolinium;  carcinoma;  kinetic  modeling; 
capillary  permeability;  blood  volume 

The  last  decade  has  seen  a  rapid  development  in  the  use  of 
dynamic  contrast-enhanced  Tj-weighted  MRI  in  medicine. 
In  tandem  with  the  technological  advances  that  have  en¬ 
abled  improved  data  acquisition,  a  number  of  investigators 
have  employed  physiological  models  to  facilitate  data  in¬ 
terpretation.  While  the  use  of  these  models  has  found 
numerous  applications  (e.g.,  in  studies  of  tumor  physiol¬ 
ogy  (1)  and  myocardial  perfusion  (2)),  and  a  number  of 
groups  have  assessed  the  potential  of  model  parameters  as 
surrogate  markers  (3,4),  little  has  been  published  that  ad¬ 
dresses  the  direct  interpretation  of  these  results.  Specifi¬ 
cally,  how  do  the  estimates  obtained  using  the  various 
models  compare  with  the  physiological  parameters  they 
purport  to  measure? 

This  is  not  a  simple  question  to  address  since  it  is  often 
difficult  to  identify  a  reliable  “gold  standard.”  Many  in¬ 
vestigators  compare  their  results  with  those  obtained  with 
positron  emission  tomography  (PET).  However,  PET 
shares  many  of  the  basic  models  employed  in  MRI  (5). 
Similarly,  data  simulation  exercises  using  Monte-Carlo 
techniques  designed  to  assess  accuracy  and  precision  in 
parameter  estimation  often  utilize  the  same  model  to  both 
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generate  and  analyze  the  data  (6,7).  In  this  way,  the  sensi¬ 
tivity  of  the  estimates  to  experimental  variables,  such  as 
noise  and  sampling  frequency,  is  assessed  but  little  is 
revealed  about  the  physiological  significance  of  the  result¬ 
ant  parameter  estimates. 

A  physiological  model  incorporating  multiple  parallel 
pathways  and  heterogeneous  flow  was  used  to  simulate 
data  of  a  realistic  nature  to  which  simplified  models  were 
fitted.  The  experiment  was  designed  to  assess  the  accuracy 
of  the  models  themselves,  not  the  quality  of  the  data  to 
which  they  are  fitted  (in  terms  of  noise,  sampling  fre¬ 
quency,  etc.),  since  data  quality  is  essentially  an  experi¬ 
mental  variable.  Furthermore,  this  study  was  restricted  to 
those  models  dealing  with  a  contrast  agent  that  diffuses 
out  of  the  vascular  space  (thereby  incorporating  capillary 
permeability  as  a  model  parameter).  Many  of  the  issues 
associated  with  the  analysis  of  data  from  susceptibility 
contrast  studies,  which  typically  assume  that  the  contrast 
agent  remains  intravascular,  have  been  investigated  in  pre¬ 
vious  studies  (e.g.,  Refs.  7  and  8).  Data  were  simulated 
using  the  distributed  pathway  model  of  tracer  kinetics 
called  the  multiple  indicator,  multiple  path,  indicator  di¬ 
lution  4  region  model  (MMID4),  made  available  by  the 
National  Simulation  Resource  (Department  of  Bioengi¬ 
neering,  University  of  Washington).  These  data  were  sub¬ 
sequently  analyzed  using  three  different  tracer  kinetic 
models,  and  the  parameter  estimates  obtained  were  com¬ 
pared  with  the  physiological  variables  used  in  the  MMID4 
model. 

Out  of  necessity,  only  a  few  models  and  the  underlying 
assumptions  that  support  these  models  can  be  examined. 
The  use  of  a  fixed  input  function  (originally  proposed  by 
Tofts  and  Kermode  (9))  is  known  to  introduce  considerable 
error  into  the  estimate  of  the  capillary  permeability-sur¬ 
face  area  product  (10).  Given  this,  the  following  analyses 
assume  that  the  arterial  input  function  has  been  measured. 
Although  techniques  to  achieve  this  are  not  yet  in  common 
use,  an  increasing  number  of  investigators  are  developing 
practical  methods  for  doing  so  (11,12). 

METHODS 

Data  Simulation 

Tissue  residue  curves  were  simulated  using  the  MMID4. 
This  model  (described  in  detail  in  Refs.  13  and  14)  ac¬ 
counts  for  flow  dispersion  and  heterogeneity,  and  includes 
capillaries  modeled  as  axially  distributed  blood-tissue  ex¬ 
change  units.  A  plasma  concentration-time  curve  was  sim¬ 
ulated  as  an  input  to  the  model.  The  amplitude  and  shape 
of  the  curve  were  based  upon  measurements  made  in  the 
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descending  aorta  of  volunteers  following  bolus  injection  of 
Gd-DTPA  (11).  Reference  data  for  comparison  with  the 
simulated  residue  curves  were  obtained  from  MR  studies 
of  human  brain,  breast,  and  prostate  tumors  (15). 

Numerous  parameters  determine  the  behavior  of  the 
MMID4  model.  Delay  and  dispersion  of  the  arterial  input 
were  set  to  zero.  Contrast  agent  in  nonexchanging  vessels 
(arteries,  arterioles,  venules,  and  veins)  may  contribute  to 
the  tissue  residue  function.  However,  such  vessels  were 
excluded  in  this  study,  as  the  tumors  of  interest  are  likely 
to  contain  a  predominance  of  newly  formed,  leaky  vessels 
(16).  Flow  heterogeneity  was  modeled  by  apportioning 
flow  through  20  parallel  pathways  according  to  a  right- 
skewed  lagged  normal  density  distribution  (14).  For  the 
sake  of  clarity  a  tissue  density  of  1  g/ml  was  assumed  for 
all  model  analyses  (17).  Other  than  mean  plasma  flow  (Fp) 
the  only  parameters  that  were  adjusted  in  the  simulations 
were  those  associated  with  blood-tissue  exchange:  capil¬ 
lary  plasma  volume  (Vp),  capillary  permeability-surface 
area  product  (PS),  and  interstitial  volume  (Ve). 

Model  Fitting 

Three  models  were  chosen  for  assessment.  Model  1,  the 
modified  Kety  model  (18),  was  essentially  applied  to  MRI 
data  by  both  Larsson  et  al.  (19)  and  Tofts  and  Kermode  (9): 

C,M  -  EFp  I  Cp(u)ex (t  -  u)Jdu  [1] 

J  0 

where  Ctjs  and  Cp  are  the  concentrations  of  contrast  agent 
in  the  tissue  of  interest  and  plasma,  respectively,  and  the 
extraction  fraction,  E  —  1  -  exp(-PS/Fp).  In  accord  with 
recent  standardization,  it  is  important  to  stress  that  the 
product  EFp  may  also  be  referred  to  as  Ktrans  (20).  Model 
2  resembles  the  first,  but  includes  a  vascular  term  (17)  and 
has  been  used  to  analyze  MR  data  in  a  number  of  studies 
(2,6): 

CHs(t)  =  VpCp(t)  +  EFp  I '  Cp(u)exp(^£  (t  -  u)Jdu.  [2] 

Model  3,  recently  described  by  St.  Lawrence  and  Lee  (21), 
potentially  enables  the  estimation  of  Fp  and  PS  separately: 

CUs(t)  =  Fp  Cp(t  -  u)du  +  EFp  Cp(u) 

^  0  —  T 

(~ef  \ 

— y—  {t  -  u  -  T)jdu  [3] 

where  r  is  the  mean  capillary  transit  time  (=  Vp/Fp ).  Model 
1  has  two  unknown  parameters  that  may  be  estimated  by 
curve  fitting:  the  product  EFp  and  the  interstitial  volume 
Ve.  Model  2  has  three  unknown  parameters:  the  product 
EFp,  the  interstitial  volume  Ve,  and  the  plasma  volume  Vp. 
Model  3  contains  four  unknown  parameters  that  may  be 
estimated  by  curve  fitting:  Fp,  F,  Ve,  and  t.  From  these,  the 
parameters  PS  and  Vp  may  be  evaluated. 


b  time  /min 

FIG.  1 .  Representative  signal-time  curves  used  for  simulations.  Tis¬ 
sue  signal-time  curves  (bold  lines)  simulated  using  an  assumed 
arterial  input  function  (faint  lines),  MMID4,  and  the  physiological 
parameters:  (a)  Fpf  0.57  ml/g/min;  Vp ,  0.06  ml/g;  PS,  0.33  ml/g/min; 
Vg,  0.45  ml/g;  and  (b)  Fpi  1.2  ml/g/min;  Vpf  0.08  ml/g;  PS, 
0.34  ml/g/min;  Ve,  0.40  ml/g.  Reference  data  (open  circles)  were 
obtained  from  studies  of  (a)  a  breast  tumor  and  (b)  a  meningioma. 

Simulated  data  were  generated  in  two  series.  The  first 
series  was  designed  to  be  representative  of  data  acquired 
from  a  breast  tumor,  and  the  second  was  representative  of 
data  acquired  from  a  meningioma  (Fig.  1).  Baseline  values 
selected  for  Vp  and  Vc  were  based  on  estimates  from  the 
literature.  These  values  were  then  held  fixed  while 
MMID4 ,  with  only  Fp  and  PS  as  free  parameters,  was  fitted 
to  the  reference  data.  Best-fit  estimates  of  Fp  and  PS  com¬ 
bined  with  the  choices  of  Vp  and  Ve  were  then  used  as 
baseline  values  for  the  simulated  data.  In  each  case  points 
were  calculated  with  a  sampling  frequency  of  one  acqui¬ 
sition  per  second  over  a  total  sampling  period  of  5  min. 
Twenty-six  tissue  residue  curves  were  produced  using  the 
range  of  parameter  values  outlined  in  Table  1  (Fig.  2).  In 
each  case  one  parameter  value  (Fp,  Vp,  or  PS)  was  modified 
from  the  baseline  values  per  curve.  Four  values  of  each 
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Table  1 

Values  of  the  Adjustable  Kinetic  Parameters  Describing  the  Behavior  of  MMID4  Used  to  Simulate  26  Different  Tissue  Residue  Curves 


Baseline 

Exps.  2-5 

Exps.  6-9 

Exps.  10-13 

(a)  Data  representative  of  a  breast  tumor 

Fp  (ml/min/g) 

0.57 

0.17,  0.37, 

0.77,  0.97 

0.57 

0.57 

Vp  (ml/g) 

0.06 

0.06 

0.0001,  0.03, 

0.09,  0.12 

0.06 

PS  (ml/min/g) 

0.33 

0.33 

0.33 

0.01,  0.17, 

0.49,  0.65 

ve  (ml/g) 

0.45 

0.45 

0.45 

0.45 

(b)  Data  representative  of  a  meningioma 

Fp  (ml/min/g) 

1.2 

0.4,  0.8, 

1.6,  2.0 

1.2 

1.2 

Vp  (ml/g) 

0.08 

0.08 

0.0001,  0.04, 

0.12,  0.16 

0.08 

PS  (ml/min/g) 

0.34 

0.34 

0.34 

0.0,  0.17, 

0.51,  0.68 

Va  (ml/g) 

0.4 

0.4 

0.4 

0.4 

The  remaining  parameters  associated  with  MMID4  were  held  fixed  (see  text).  A  representative  selection  of  the  resultant  tissue  residue 
curves  is  shown  in  Fig.  2. 


parameter  were  chosen  in  addition  to  the  baseline  values; 
thus,  13  curves  were  generated  per  series.  Each  of  the  three 
models  in  turn  was  fitted  to  the  26  simulated  data  sets  and 
estimates  of  the  parameters  were  compared  with  the  val¬ 
ues  used  in  MMID4  to  generate  the  residue  curves. 

RESULTS 

Data  Simulation 

The  simulated  data  exhibit  many  of  the  characteristics  of 
the  experimental  data  (Fig.  1).  Indeed  the  simulated  data 
match  the  example  data  very  well  despite  the  arterial  input 
function  used  (a  single  simulated  input  not  determined  in 
individual  subjects).  The  physiological  parameters  influ¬ 
ence  the  residue  curves  in  different  ways  (Fig.  2).  The 


initial  slope  of  the  residue  curve  is  determined  principally 
by  the  vascular  parameters.  The  early  curvature  and  peak 
observed  in  the  data  is  strongly  influenced  by  PS.  There  is, 
however,  considerable  overlap  of  these  effects. 

Model  Fitting 

Models  1  and  2  fitted  the  simulated  data  in  a  reproducible 
manner.  Parameter  estimates  obtained  from  these  fits  are 
displayed  in  Fig.  3.  Parameter  estimates  obtained  using 
model  3  were  not  always  unique,  and  often  depended  on 
the  starting  guesses  used  in  the  fitting  routine.  This  im¬ 
plies  a  significant  correlation  between  parameters  and  the 
existence  of  local  minima  (22).  To  improve  the  stability  of 
the  regression  procedure,  an  approach  similar  to  that  de- 


FIG.  2.  Influence  of  the  physiological  parameters.  The  shape  of  the  signal-time  curve  is  modified  by  changes  in  the  values  of  individual 
parameters.  The  baseline  curve  (bold  line)  is  representative  of  data  acquired  from  a  breast  tumor.  The  influence  of  changes  in  (a)  plasma 
flow  (Fp),  (b)  permeability-surface  area  product  (PS),  and  (c)  plasma  volume  (Vp)  are  shown.  The  parameter  values  used  to  generate  the 
curves  are  detailed  in  Table  la. 
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FIG.  3.  Parameter  estimates  as  a  function  of  the  “true”  values.  Estimates  of  the  parameters  (a)  EFp)  (b)  Vp,  and  (c)  Fp  and  PS  obtained  by 
fitting  model  1  (squares),  model  2  (circles),  and  model  3  (triangles)  to  the  simulated  data  (Table  la).  The  open  triangles  in  c  represent 
estimates  of  Fp,  the  filled  triangles  represent  estimates  of  PS,  and  the  dashed  lines  in  all  figures  represent  lines  of  identity.  Similar  plots  were 
obtained  using  the  parameters  of  Table  1b.  a:  Note  the  failure  of  model  1  at  low  EFP  (0.01  ml  min“1g“1).  The  model  1  estimate  is  off-scale 
at  0.82  ml  min_1g_1. 


scribed  by  Henderson  et  al.  (22)  was  adopted.  Multiple  fits 
were  performed  using  a  range  of  starting  points  (i.e.,  dif¬ 
ferent  initial  estimates  of  t).  The  final  solution  selected 
was  that  which  minimized  the  sum  of  squared  differences 
between  the  data  and  the  model  fit.  Representative  results 
are  displayed  in  Fig.  3.  Though  the  results  in  general  were 
convoluted,  certain  trends  were  apparent.  In  particular,  Vp 
was  consistently  underestimated  by  both  model  2  and 
model  3  (by  between  2%  and  96%).  The  parameters  of 
model  3  became  less  accurate  at  low  flow  or  high  perme¬ 
ability,  where  PS  tended  to  be  overestimated  (to  an  ex¬ 
treme  of  152%)  and  Fp  was  slightly  underestimated  (up  to 
34%). 

DISCUSSION 

Increasing  numbers  of  investigators  are  employing  models 
to  analyze  their  dynamic  contrast-enhanced  data  (20).  It  is 
important  that  the  techniques  used  are  fully  understood  so 
that  they  may  benefit  other  investigators.  Part  of  this  pro¬ 
cess  is  to  investigate  the  physiological  interpretation  of 
these  parameters.  Typically,  this  would  be  achieved  by 
comparison  of  the  MRI  results  with  a  “gold  standard.” 
However,  suitable  techniques  for  comparison  are  difficult 
to  find.  PET  has  often  been  quoted  as  a  gold  standard 
technique  for  measurement  of  blood  flow  (8),  but  the  tech¬ 
niques  used  to  interpret  PET  data  are  very  similar  to  those 
now  used  in  MRI  and  can  suffer  from  similar  pitfalls  (5).  In 
this  study,  a  realistic,  established  model  has  been  used  to 
simulate  tissue  residue  curves  that  closely  match  existing 
experimental  data.  In  this  way,  data  were  generated  with 
known  physiological  variables  that  were  independent  of 
the  models  subsequently  used  to  interpret  them. 

Data  Simulation 

MMID4  is  clearly  capable  of  generating  residue  curves  that 
closely  match  experimental  data  acquired  by  MRI  (8,14,23) 
(Fig.  1).  This  should  not  be  interpreted  as  suggesting  that 


MMID4  fully  describes  the  physiology  of  the  tumors  stud¬ 
ied.  It  is  simply  a  step  towards  that  goal.  However,  MMID4 
does  incorporate  parameters  of  particular  interest:  Fp,  PS, 
Ve,  and  Vp,  and  each  influences  the  residue  curve  over 
differing  temporal  ranges  and  in  subtly  different  ways  (Fig. 
2).  Fp ,  for  example,  determines  the  initial  slope  of  the 
residue  curves  simulated  in  this  study.  PS  has  a  more 
convoluted  impact  depending  upon  the  ratio  PS :Fp.  Cer¬ 
tainly,  the  curvature  of  the  residue  curve  about  the  initial 
peak  (Fig.  2)  is  strongly  influenced  by  PS.  Ve  has  perhaps 
the  clearest  affect  upon  the  residue  curve.  Increasing  Ve 
tends  to  delay  the  rise  of  the  curve  to  its  equilibrium 
amplitude.  Finally,  the  impact  of  Vp  is  seen  throughout  the 
time  course.  It  is  the  overlap  in  influence  of  each  variable 
that  complicates  the  issue  of  parameter  estimation. 

Model  Fitting 

A  model,  by  definition,  is  a  practical  simplification  of  the 
true  tissue  physiology;  as  a  result  it  is,  to  some  extent, 
flawed.  It  is  therefore  not  surprising  that  the  simple  mod¬ 
els  tested  were  unable  to  provide  accurate  estimates  of  the 
physiological  parameters  in  all  cases.  Model  3  bears  the 
closest  resemblance  to  MMID4  but  clearly  fails  to  describe 
the  more  complex  model  in  certain  circumstances.  In  par¬ 
ticular,  model  3  systematically  underestimates  Vp  and 
overestimates  PS  under  conditions  of  low  flow  or  high 
permeability  (Fig.  3).  This  may  be  a  consequence  of  flow 
heterogeneity  (14),  an  effect  all  three  models  overlook. 
Flow  heterogeneity  differs  between  regions  and  individu¬ 
als  and  may  not  remain  constant,  as  assumed  in  this  study, 
but  could  change  as  a  function  of  pathology  (8).  The  effect 
of  additional  disrupting  influences  on  the  input  function — 
delay  and  dispersion — have,  for  the  sake  of  conciseness, 
not  been  considered  in  this  study  since  their  influence 
upon  parameter  estimation  has  been  studied  elsewhere  (7). 

If  the  direct  vascular  contribution  to  the  signal  is  ignored 
(model  1),  EFp  (also  known  as  K lrans  (20))  is  overestimated 
(by  up  to  54%  (Fig.  3a)).  Indeed,  under  certain  circum- 
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stances  model  1  breaks  down  completely.  However,  to 
measure  the  vascular  signal  contribution  it  is  essential  to 
measure  the  vascular  input  and  acquire  data  with  a  high 
temporal  resolution  (6).  The  choice  of  model  3  for  data 
analysis  requires  a  sampling  interval  less  than  the  mean 
transit  time  of  the  tracer.  Once  the  sampling  interval 
reaches  the  mean  transit  time,  model  3  may  be  simplified 
to  a  form  similar  to  model  2  (21).  It  is  clear  from  Fig.  3  that 
the  use  of  model  2,  despite  the  inclusion  of  a  vascular 
contribution,  leads  to  consistent  underestimation  of  Vp 
and  a  slight  overestimation  of  EFp.  These  findings  mirror 
the  results  of  a  previous  study  (21).  Interestingly,  the  ratio 
EFp:Vc  may  be  estimated  with  an  error  of  less  than  10%, 
which  is  consistently  less  than  the  error  in  the  estimate  of 
EFp  alone  (data  not  shown).  This  ratio,  otherwise  referred 
to  as  kep  (20),  may  be  estimated  without  a  tissue  7\  mea¬ 
surement  (19),  and  has  been  identified  as  a  useful  param¬ 
eter  for  differential  diagnosis  in  the  study  of  primary  breast 
tumors  (24). 

Another  important  issue,  which  is  not  always  apparent 
to  the  investigator,  is  the  uniqueness  of  the  model  solution 
(23).  Many  of  the  estimates  obtained,  using  model  3,  varied 
as  a  function  of  the  starting  values  used  for  curve  fitting. 
The  array  of  local  minima  result  from  parameter  correla¬ 
tions,  i.e.,  different  combinations  of  the  parameters  pro¬ 
duce  very  similar  solutions.  Parameter  correlations  may  be 
examined  using  sensitivity  functions  (13)  or  by  examina¬ 
tion  of  the  fit  correlation  matrix.  Parameter  estimation  is 
simplified  if  the  complexity  of  the  fit  can  be  reduced.  It  is 
possible  to  use  MMID4  for  data  analysis  (14,23),  but  the 
issue  of  parameter  uniqueness  is  further  complicated  by 
the  surfeit  of  possible  parameter  combinations.  Jerosch- 
Herold  et  al.  (23),  studying  myocardial  blood  flow  and 
capillary  permeability  in  a  porcine  model,  found  that  they 
were  unable  to  obtain  unique  estimates  of  blood  flow, 
volume,  and  PS  using  an  extracellular  contrast  agent  (Gd- 
DTPA).  They  were  able  to  assess  blood  flow  using  an 
experimental  intravascular  tracer,  as  all  the  extravascular 
parameters  were  eliminated  as  degrees  of  freedom  of  the 
model.  Even  under  these  circumstances,  it  was  necessary 
to  fix  the  value  of  numerous  additional  MMID4  parameters 
at  literature  values  (23).  It  is  likely  that  future  studies  may 
benefit  from  the  combined  use  of  intravascular  and  extra¬ 
cellular  tracers  in  multiple-indicator-type  modeling  (13). 
Indeed,  model  3  has  been  tested  with  two  different  con¬ 
trast  agents  (22)  in  a  canine  model  of  breast  cancer.  Fur¬ 
thermore,  Henderson  et  al.  (22)  recognized  the  confound¬ 
ing  impact  of  parameter  correlations  upon  parameter  esti¬ 
mation.  To  overcome  this  problem,  curve  fitting  was 
performed  in  a  serial  fashion  using  a  succession  of  starting 
points  for  the  fit  (22),  a  policy  adopted  in  this  study. 

Implications  for  Tracer  Kinetic  Studies 

It  is  important  at  the  outset  of  a  study  to  identify  the  most 
relevant  parameters  to  measure.  Without  an  arterial  input 
function  and  with  limited  temporal  resolution,  the  Tofts 
and  Kermode  (9)  and  Larsson  et  al.  (19)  models  have 
proved  very  useful.  In  tissues  with  a  negligible  blood  vol¬ 
ume  they  may  also  provide  good  estimates  of  EFp.  How¬ 
ever,  with  a  measured  arterial  input  function  and  rapid 
data  collection,  model  2  may  be  a  more  appropriate  model 


FIG.  4.  Measurement  of  the  arterial  and  tissue  concentrations  in 
vivo.  Given  an  appropriate  imaging  sequence  the  arterial  input  func¬ 
tion  and  the  tissue  residue  curve  can  be  measured  in  vivo.  These 
data  were  acquired  from  regions  of  interest  encompassing  the 
sagittal  sinus  (open  squares)  and  a  meningioma  (open  circles). 
Calculation  of  the  concentration  of  Gd-DTPA  (mM)  was  based  upon 
precontrast  T-,  measurement  (12). 


for  clinical  data  evaluation.  In  either  case  the  investigator 
must  appreciate  the  assumptions  and  implications  of  us¬ 
ing  these  models.  In  particular,  neither  model  provides 
separate  estimates  of  perfusion  and  permeability.  These 
parameters  are  inextricably  combined  in  the  EFp  product. 
The  model  (model  3)  of  St.  Lawrence  and  Lee  (21)  offers 
the  possibility  of  assessing  blood  flow  and  extraction  si¬ 
multaneously,  but  may  require  carefully  constrained  su¬ 
pervised  data  fitting  (22),  additional  measurements,  or 
multiple  tracers  to  overcome  problems  with  parameter 
correlations  (23).  In  all  cases  it  is  essential  to  collect  ap¬ 
propriate  data  of  the  highest  quality  if  accuracy  and  pre¬ 
cision  are  to  be  attained. 

As  a  practical  example  of  these  choices,  the  models  were 
each  fitted  to  recent  experimental  data  (input  and  tissue 
residue  curves)  acquired  in  a  meningioma  (12),  shown  in 
Fig.  4.  These  data  were  simultaneously  acquired  every 
5.1  s  in  the  sagittal  sinus  and  tumor  using  a  ^-weighted 
3D  gradient-echo  sequence  (TR  =  4.3  ms,  TE  =  1.1  ms,  flip 
angle  =  35°)  on  a  1.5  T  MR  system.  The  signal  from  the 
sagittal  sinus  was  used  as  input  to  the  models  (a  vein  was 
chosen  for  simplicity  and  is  merely  illustrative)  and  the 
effective  sampling  interval  was  halved  by  linear  interpo¬ 
lation.  The  following  parameter  estimates  were  obtained 
using  model  3:  blood  (rather  than  plasma)  flow,  Fb  = 
0.25  ml  min-1g-1;  blood  volume,  Vh  =  0.02  ml  g_1;  Ve  ~ 
0.43  ml  g-1;  and  PS  =  0.14  ml  min-1g-1.  It  is  unclear 
whether  the  interpolation  step  was  required  since  the  cap¬ 
illary  transit  time,  t,  was  estimated  to  be  4.8  s.  Based  on 
the  results  of  the  previous  simulation  exercise,  PS  and  Fb 
may  be  somewhat  overestimated  and  underestimated,  re¬ 
spectively.  Furthermore,  the  estimate  of  Vb  is  likely  to 
represent  a  significant  underestimate.  However,  the  prod¬ 
uct  EFb  (=  0.16  ml  min-1g-1)  is  likely  to  be  accurately 
estimated,  and  these  results  compare  well  with  estimates 
obtained  using  model  2:  EFb  =  0.16  ml  min-1g-1;  Vb  = 
0.02  ml  g"*1;  and  Ve  =  0.44  ml  g"1.  Predictably,  the  estimate 
of  the  EFb  product  made  by  model  1  was  larger  (0.18  ml 
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min^g"1}  but  a  comparable  value  was  obtained  for  Ve 
(0.44  ml  g“a).  Unless  an  independent  estimate  of  blood 
flow  or  permeability  is  essential,  model  2  may  be  an  ap¬ 
propriate  choice  for  analysis  of  these  data  since  it  requires 
neither  data  interpolation  nor  supervised  fitting. 

CONCLUSIONS 

The  use  of  models  for  tracer  kinetic  analysis  potentially 
enables  the  investigator  to  interpret  MR  data  in  physiolog¬ 
ical  terms.  However,  both  the  measurement  process  and 
the  data  analysis  introduce  a  degree  of  uncertainty  into 
parameter  estimation.  Simplified  models,  some  of  which 
are  in  common  use,  may  not  provide  accurate  estimates  of 
the  physiological  parameters  of  interest.  It  is  important  to 
address  these  uncertainties  because  it  is  only  after  the 
parameters  of  the  model  have  been  accurately  estimated 
that  the  diagnostic  or  prognostic  efficacy  of  specific  quan¬ 
titative  physiological  parameters  can  be  properly  tested. 
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